On the effectiveness of threshold resummation away from hadronic endpoint 
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We parameterize the enhancement of threshold effects away from hadronic endpoint that arise 
due to the steeply falling nature of parton distribution functions, within the context of soft-collinear 
effective theory. This is accomplished in a process-independent way by directly linking the charac- 
teristic scale of soft and collinear radiation, A, to the shape of the pdfs. This allows us quantify 
the power corrections to partonic threshold resummation as a function of the invariant mass and 
rapidity of the final state. In the context of SCET, being able to compute A in a process-independent 
manner allows us to determine the correct scale for threshold resummation after integration with 
the pdfs, without any additional procedure. 



INTRODUCTION 



There has been much effort to increase the accuracy and precision of predictions for obscrvablcs at hadron colliders. 
The standard technique to improve the accuracy of calculations is to add fixed orders in perturbation theory. Many 
observables have been calculated to next-to-leading order (NLO) in QCD, and for some next-to-next-to-leading order 
(NNLO) has been achieved [H-fllj. Frequently, observables contain large ratios of scales, usually due to experimental 
cuts or the presence of several mass scales in the process. In such cases, one can increase the precision and accuracy of 
theoretical calculations by rcsumming the large logarithms of these ratios to all orders in perturbation theory. Standard 
techniques allow for the resummation of these logarithms at next-to-leading logarithmic order (NLL) and beyond (]~2| — 
Qh [HI, while recent advances in effective th eory methods fl5l - [l8j provide a systematically improvable alternative to 
. resummation in hard-scattering processes fl9j . 

In this paper, we will focus on so-called threshold logarithms. It is well known that, beyond leading order, partonic 



cross sections contain terms of the form 

(N 



where m < 2n — 1 and z = s m j n /s is a measure of excess radiation in the process. These terms are due to collinear 
. ■ and soft singularities in the real emission diagrams of the perturbative series and can be resummed with the same 
techniques used for other large logarithms. Threshold resummation has been applied to several Standard Model 
processes, including (inclusive) Drell-Yan [H, [H, HU , prompt photon [22|, [23|, Higgs production [2^ - [28| . and dijet 
and heavy-particle production (29j - [34| . 

It is important to understand when threshold logarithms are phenomenologically relevant. The partonic endpoint 
(where logarithms of z become large) cannot be observed, since one integrates over the partonic variable z when the 
partonic cross section is convolved with parton distribution functions (pdfs). The exception is near the hadronic 
endpoint, where the partonic and hadronic endpoints are equivalent. As an example, consider Drell-Yan, where the 
' rescaled leptonic center of mass energy is given by r = m 2 u /s. The hadronic cross section is then 



-j^ = J dzC(T/z)a(T,z) , 



where 



flog f 

C(z) = I dY 

log 7 



h(x iy) h (f z e-r) +h (^ e -r ]h(xl ^ 



(2) 



(3) 



is the parton luminosity. From this expression, one can immediately sec that a measurement of Drell-Yan close to the 
hadronic endpoint (r — > 1) also forces the partonic center of mass energy z to its endpoint, such that these logarithms 
become important. However, away from the hadronic endpoint, the integration variable z is not forced to one, such 
that it is not clear that the threshold terms given in Eq. ([1} should dominate over terms which are, for example, 
polynomial in z. 

By comparing resummed cross sections with available fixed order calculations it has been noted that, in certain 
cases, the threshold terms at 0(a s ) give the dominant effect of the full NLO correction, even far away from the 
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hadronic cndpoint [U [22|, [H| . This observation has been used to suggest that threshold resummation is effective 
not only in the hadronic endpoint, where it is clearly necessary, but also for much lower values of r. One possible 
explanation for this is the shape of parton luminosities [13,|35|. For steeply falling parton luminosities, the integration 
in Eq. ([5]) is dominated by the smallest values of t/z, which is the region z — > 1. By assuming a simple, analytical 
form for the parton luminosity, one can make this argument more precise (2lj . However, the corrections to threshold 
resummation applied away from hadronic endpoint have not yet been studied without assuming a model for parton 
luminosities. 

The goal of this paper is to provide a model- independent, quantitative measure of the corrections to partonic 
threshold resummation. Using only the shape of the pdfs, we define a parameter A that can be used as the expansion 
parameter in SCET. This implies that corrections to the threshold resummation are power suppressed in A. Since 
we will use the measured pdfs to calculate A, our definition does not require the assumption of a particular form 
for the parton luminosities. Thus, A can be viewed as a mo del- independent definition of the steepness of a parton 
luminosity. We will determine the numerical size of A for various parton luminosities as a function of the hadronic 
threshold variables r and Y. In addition to quantifying power corrections, A allows us to determine the scales relevant 
to threshold resummation in SCET after convolving the partonic cross section with the pdfs. This avoids integrating 
over unphysical regions (a consequence of setting the scales before integration) or introducing a process-specific 
procedure to calculate the scales. 

This paper is organized as follows. In Section [3J we will define the relevant expansion parameter in threshold 
resummation and show that it has the expected behavior for a simple choice of the pdfs. In Section [31 we will 
determine the numerical values of this expansion parameter for the actual pdfs, and present our conclusions in 
Sectional In Appendix [3] we show the consistency of the factorization theorem assuming only small values of A, but 
keeping an arbitrary functional form for the pdfs, while in Appendix [B] we discuss the scaling in A of integrals against 
the pdfs. 



2. DEFINITION OF STEEPNESS 



The only assumptions one has to make to allow for the resummation of threshold logarithms is that 1 — z <C 1, 
where z is appropriately defined, along with the statement that pdfs factorize from the partonic cross section [36l.l37|. 
We will use the definitions 



z = — 



(4) 



with 



N 



Q + Y,PJ) ■ ( 5 ) 

i 

Here q denotes the sum of the momenta of all non-strongly interacting particles, and j»j is the momentum of the ith 
jet. This momentum is defined in terms of pj and r/j of the jet as 

P.J = (Pj cosh 77,7, pj,p T j sinhrjj) . (6) 

Thus, the jet 4-momentum is reconstructed from the measured transverse momentum and rapidity, assuming zero 
mass. 

A generic hadronic cross section can be written as 



da = J ^ dYj 2Y dzf 1 ^- z e Y Jf 2 ^- z e- Y )da(z,T,Y) + (l^2). (7) 

One can use soft-collinear effective theory (SCET) to analyze the partonic cross section, with corrections being 
suppressed by powers of 1 — z 

da(z, r, Y) = da(z, r, r) SCE T + 0(1 - z) . (8) 

As was shown in [38j . the SCET partonic cross section factorizes, and renormalization group (RG) equations can be 
used to resum threshold logarithms. This work made no assumptions about the kinematics of the event other than 
1 — z ~ A 2 , where A is the usual SCET power counting parameter, which relates the hard, collinear and soft scales 
via ii c jii h ~ Us/ He ~ A. 



3 



Prom Eq. (JT]) it is clear that z is forced to one if the kinematics force us close to hadronic threshold, r 1. In this 
case, 1 — z -C 1, which means that the power corrections in Eq. ((5]) are small. Away from hadronic threshold, however, 
the partonic cross section gets integrated over values of z for which 1 — z ~ O(l). For generic parton luminosities, 
threshold resummation will have 0(1) power corrections when applied away from hadronic threshold, if the region 
1 — z ~ 0(1) is large. 

If the parton luminosities are dominated by the region where 1 — z <C 1, however, the SCET expression gives 
the dominant contribution to the hadronic cross section. We can quantify the corrections due to this by defining a 
parameter 



e(A,T,Y) 



f re -2\Y\ dz 






sLx» dz 1 


A(v^)l 


h(^~ Y )\ 



(9) 



In other words, the parameter e measures how much of the parton luminosity is contained in the region re _2 ' y ' < 
z < 1 — A 2 compared to the region 1 > z > 1 — A 2 . This parameter obviously depends not only on our choice of A, 
but also on the hadronic center of mass energy r and the total rapidity of the event Y. Note that for fixed values of 
r and Y, e is as a monotonically decreasing function of A. 
Using this definition, we can write 



dcr 



dY 



-illl 



dzh 



1-A 2 



da(z,T,Y) SC ET + 0(e,A 2 ), 



(10) 



where the 0(A 2 ) correction comes from the fact that 1 — z in Eq. (jlpp is bounded by A 2 . The total correction to 
threshold resummation can therefore be estimated as max(e, A 2 ). However, since e decreases with increasing A 2 , the 
total correction is minimized if we choose e = A 2 . The corrections are then 0(A 2 ). This also allows us to define A 2 
through the relation 



A 2 = 



j Te -2\Y\ dz 




|/ a (vT«" y )l 


sL\' dz 1 







(11) 



where A 2 = A 2 (r, Y). We can also define A 2 for a single pdf 



A 2 



t ^ ^ |/ (f ) 
Ji-\? dz |/(f)l 



(12) 



where here A 2 = A 2 (x). This definition will be useful for determining if endpoint Altarclli-Parisi functions are valid 
when evolving the pdfs. 
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FIG. 1: Left: The parameter e, defined in Eq. ([9]), illustrated as the ratio of two areas. Right: For power corrections that go as 
max(e, A 2 ), there is a choice of A 2 = Aj| nal such that the corrections are minimized. This choice, called A 2 in the text, is defined 
in Eqs. {TTJ and (fl2]t. 



We illustrate our definitions in Fig. [TJ with a simple example function f(x). The graphic on the left shows e as the 
ratio of two areas: the region D, where 1 — z < A 2 and the region N, where 1 — z > A 2 . For N <C D and A 2 <C 1, we 
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see that the area under the curve is dominated by a region where 1 — z <C 1. The right graphic shows that Eq. (jT2|) 
has a solution where power corrections are minimized, labelled A^ nal . 

Independent of the form of the pdfs, we can see that A| nal goes to as r (or x) goes to 1, as expected. First, note 
that t — y 1 forces Y — > 0, which sets the lower bound of the integral in the numerator of Eq. © to r. Next, we can 
examine the limits on A 2 . As A 2 goes to 0, e goes to oo (here we've assumed that the pdf implicitly contains a 
function that sets it to if the argument is greater than 1 or less than 0). Similarly, as A 2 goes to 1 — r, e goes to 
0. Therefore, the crossover point, where e = A 2 , is between and 1 — r, which means A| nal — > as r — > 1. For the 
remainder of the paper, we will omit the subscript final on A 2 . 

Note that these definitions do not assume anything about the form of the pdfs and are therefore completely model 
independent. Given that pdfs are not observable quantities, many pdf sets do not constrain them to be positive. 
For this reason, we have used absolute values of the pdfs in our definition. Of course, the value of A 2 depends on 
the functional form of the pdfs, but this can be determined numerically using any of the available pdf sets. We will 
present the numerical results for the power counting parameter A 2 in Section [3] 

We end this section by showing the behavior of A 2 for a particularly simple form of the pdfs. This will illustrate our 
proposal with an analytic example and should convince the reader that corrections to threshold rcsummation behave 
as expected. The assumed form of the pdfs docs not represent their observed shape and is only chosen so that the 
analytic results can be easily studied. All numerical results presented in Section [3] will use the measured pdfs. 

Consider the simple model 

fi{x)=x—, (13) 
where the pdf becomes steeper as a increases. Given this form, Eq. becomes 

2 (1 - A 2 ) 1+ ( Ql + a2 )/ 2 - ( Te - 2 l^l) 1 + ( a i+ a 2)/2 
A = 1 - (1 - \2y + ( ai +a 2 )/2 ' ( 14 ) 



To simplify this equation, we define 



r 



= AX 2 , with A=l + ^^. (15) 



The equation can then be written as 

We can clearly see that r — > in the hadronic endpoint r — > 1 (for which Y — > 0), independent of the value of the 
parameter A. However, for large values of A, and therefore steep pdfs, we can use the approximations 1 + r/A w 1 
and (1 - r/A) A w e~ r to find 

r = W{A) , (17) 

where VF(^4) denotes the product logarithm of A. Since the IF(yl) only grows logarithmically with A, we find the 
desired result that A 2 — s- as A —> oo. 



3. NUMERICAL RESULTS FOR A 2 



Now that we have defined a measure of steepness, it is a simple matter to determine numerically what A 2 is for the 
physical pdfs and parton luminosities. In this section, we will give our results and discuss their implications for the 
relevance of threshold rcsummation. As discussed above, the expansion parameter for threshold resummation is given 
by A 2 , and we make the assumption that this expansion becomes reliable for A 2 < 0.25. Note that while this value is 
chosen somewhat ad-hoc, it serves to assess the importance of threshold resummation. 

There are two different types of objects for which steepness is a necessary ingredient of our factorization theorem. 
The first is individual pdfs, where A 2 C 1 allows us to use endpoint AP kernels in determining the consistency of 
our factorization theorem. The second is parton luminosities, in the form of two pdfs which are convolved in Y 
and z. Since it is possible to be differential in Y up to power corrections in A 2 (see Appendix |A")) . we will focus 
on the case where Y is fixed and z is the only convolution variable. We will present plots using the CTEQ6.6 [H|, 
MSTW2008NLO Hj and NNPDF2.0 [HJ pdf sets, with a default renormalization scale of 1 TeV. To study the 
uncertainties in our predictions for A, we use the 100 replica set provided by NNPDF2.0. The central value is given 
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by the median of the results, while the error bars show the range where 68% of the points lie, with equal number of 
points on either side of the median. 

In Fig. [2] we present the value of A 2 for a single pdfs as a function of x, as defined in Eq. (fT2j) . which characterizes 
how important the non-singular terms in the AP evolution kernels are. In each plot we show the result for CTEQ, 
MSTW and NNPDF sets, with errors presented for the NNPDF set. One clearly sees that the value of A 2 is decreasing 
as x — > 1, as expected, and that A 2 is small not only in the region where 1-j;<1 but for smaller values of x as 
well. Irrespective of the type of parton, the value of A becomes smaller than 0.25 for x ~> 0.3, implying that in that 
region the singular AP evolution kernel becomes a reasonable approximation to QCD, even though 1 — x = 0(1). 
The feature in the sea-quark distributions for x ~ 0.5 is due to the steepness falling off in the measured pdfs, giving a 
rising value of A 2 , which then gets forced to zero due to x approaching 1. In the lower right plot of Fig. [21 we illustrate 
the dependence of steepness on the renormalization scale of the pdfs. Between 100 GeV and 10 TeV, there is at most 
an 0(5%) variation, while for most values of x the variation is 0(1%). In this plot, we only present the results for 
the down quark, however, the effect is similar (0(1 — 5%)) for all other pdfs. 

Next, we study the value of A 2 as obtained from the convolution of two pdfs, as defined in Eq. (fTT|) . The results are 
shown in Fig. [3] for Y = 0. While the distributions in general do not peak exactly at Y — 0, the peak is sufficiently 
close to this value for these plots to illustrate the observed behavior. Again we show the result for CTEQ, MSTW 
and NNPDF sets, with errors presented for the NNPDF set. The value of A 2 becomes less that 0.25 for r > 0.05 - 0.1. 
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FIG. 3: The value of A 2 for a selection of parton luminosities as a function of r, as defined in Eq. 



This implies that threshold resummation is not limited to the region where 1-t < 1, however, for most values 
of phenomcnological interest, the power corrections are potentially large. In the lower right hand side of Fig. [3] 
we show the variation of steepness with Y, presented as a function of Y max = —(1/2) In r. While this variation is 
considerable, one should keep in mind that the majority of the phase space is close to Y = 0. The variation of A 2 
with renormalization scale is not shown for the parton luminosities, but is similar to the result for the single pdfs. 

4. CONCLUSIONS 

In this paper we presented a model-independent, quantitative measure of the power corrections to partonic threshold 
resummation. Previously, any quantitative study of the power corrections associated with applying hadronic threshold 
resummation away from the true endpoint necessitated simplified assumptions about the functional form of the pdfs. 
The expansion parameter of threshold resummation is given by a parameter A 2 , which can be defined unambiguously 
through integrals over pdfs. In SCET, the parameter A 2 is related to the scale of soft and collinear radiation, such that 
the scales /j, c and fi s can be chosen after convolving the partonic cross section with the pdfs. This allows us to avoid 
integrating through regions where the pcrturbative expansion is ill-defined, without having to define an arbitrary 
procedure for determining scales after integration. We have shown analytically that the value of A approaches zero 
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cither in the true hadronic endpoint or in the limit of infinitely steep pdfs. Our numerical results indicate that while 
threshold resummation can be justified away from the region where 1 — r <C 1, we expect the corrections to become 
sizable for r < 0.05 - 0.1. 
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Appendix A: Consistency of threshold resummation 



As discussed in the introduction, threshold resummation can be derived using effective field theory methods. Using 
SCET, and assuming 1 — z ~ A 2 -C 1 as well as da = / <g> / <E> da, one can derive a factorization theorem for a generic 
differential cross section. We refer the reader to [H| for the details of this derivation. One check of this factorization 
theorem is consistency, which is the fi independence of the factorized result. In (38j we showed the consistency of our 
factorization theorem using the partonic definition of the pdfs = (5(1 — x). Since we have just shown that 

away from the true hadronic endpoint it is the functional form of the pdfs that defines the expansion parameter A, we 
wish to repeat this proof using the full functional form, subject only to the constraint of large steepness (or small A). 

We start from the derivative with respect to In /j, of the factorized differential cross section derived in [38| , together 
with the anomalous dimensions for the hard, jet, and soft functions, given in the same reference. Using the result of 
Appendix [B] for the Altarelli-Parisi splitting kernel, and performing trivial integrations and cancellation, we find 



da 
din// 



oc 



&Yh{^e Y )h{^e- Y ) 



T 2 In 



T 2 ln 



-4+E T ^ ln — 



rii ■ rij 



V T 2 ln^ 

«E{jcts} 




fx 



Te |21' 
1 



-v h{y/re Y )(l-w 2 )+ J Te \^\ ~~ fx {\/re Y ) / 2 («Jre Y ) (l-z)+^ 
The first line in Eq. (|Al|) can be shown to vanish using the identity 



dz 



i 



(id) 



E 



i£{partons} 



Tf In — 



(Al) 



(A2) 



The first term in the second line can be made identical to the second term through the substitution Y' =Y—\ log w\ , 
followed by the relabeling Y' — > Y. Thus these two terms cancel, and using a similar substitution Y' = Y + i logu>2, 
with the same relabeling, the third line vanishes as well. Thus, we find 

^ = + O(A 2 ), (A3) 

and therefore consistency of the factorization formula. Note that in both cases Y' = Y + 0(A 2 ), which means we can 
be differential in Y, up to power corrections in A 2 , and still have consistency. 



Appendix B: General power corrections 

It is instructive to understand the power corrections and scaling of various test functions integrated against a steep 
distribution, f(z), where steep is defined as 

Jz ""» - A < 1 . (Bl) 
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We define all scaling relative to fi_ X 2 dz f(z), which can be taken to be 0(1) without loss of generality. 

The simplest choice for a test function is a non-singular polynomial of z, g(z). In general, we can express g as a 
Taylor series about z = 1, 

where we assume all g(")(l) are 0(1). We can evaluate the integral of f(z)g(z) by dividing the integration region into 
two parts, 

dz f(z)g(z) = £ A dz f(z)g(z) + £ ^ dz f(z)g(z) . (B3) 

In the region z m - m < z < 1 — A 2 , the integral is 0(A 2 ), since /(z) and f(z)g(z) are of the same order and the integral 
of /(z) over this region is 0(X 2 ). For 1 — A 2 < z < 1, 



' dzf{z )g {z ) = J2^ 1 dzf(z)(z-lY 



ESr/ dz /( z)x(AT. (B4) 



n=0 

In this region, the integral is dominated by the zeroth order term in the power expansion. Therefore, we can write 
the integral of g{z) against the distribution /(z) as 

dzf(z)g(z)=g(l) [ dzf(z)+0(X 2 ). (B5) 

Jl-A 2 

A second possible test function is a (5-function at z = 1. The integral is trivially /(l), but we would like to know 
the scaling of this result. Assuming / is monotonically increasing from z m j n to 1, 

l 

dz/(z)</(l)xA 2 . (B6) 

l-A 2 

Since this integral is 0(1), we can say that /(l) is 0(A~ 2 ). 

The last class of test functions we will consider is log"(l — z)/(l — z) plus functions, multiplied by a non-singular 
function g{z). When integrated against /(z), this gives 

1 dzf(z)g(z)(^^ / + 

1 dziog-q-z)^-^ 1 ) - /(imp pd^a-*) 

min 1 — z 'Jo 1 - z 

1 dz log«(l - ,)/(«)*(') -/(W 1 ) + r 2 dz l 0g «(l - z) l^m + f(l)g(l)^f . (B7) 

l-A 2 



1 — z / 1 — z ' n + 1 



Since /(z) is at most 0(A 2 ) over the interval z m - m < z < 1 — A 2 , the second term is clearly power suppressed relative 
to the final term. Focusing on the first term, if we express g as a Taylor series about z = 1, we see that the zeroth 
order term is again dominant. This gives 

1 d2/( W!2£(i^A =9(1) /' cU/w f!2S^£)\ +oOT . (B8) 



•/l-A 2 V / + 



Note that Eq. (|B7j) shows that the scaling of this integral is 0(A~ 2 log n+1 A 2 ) and therefore, for small A 2 , threshold 
effects are important. 

The Altarelli-Parisi (AP) kernels, Pij(z), are polynomial functions of z for i ^ j and a combination of 5- functions 
and plus functions for i = j. This means that for steep pdfs (assuming all A 2 are small), the Pij for i ^ j arc 



9 



suppressed in A 2 compared to Pa. Moreover, any polynomial functions of z in Pa(z) can be evaluated at z = 1, up 
to power corrections in A 2 . The result is endpoint AP kernels, given at NLL accuracy by 
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